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Abstract 

The  report  outlines  the  methodology  utilized  and  results  acheived  under  the 
U.  S.  Army  Research  OfTice  grant  number  DAAL03-90-G-0221.  The  work  covered 
two  main  areas  of  resca''ch — the  aerodynamics  of  rotor  blades  including  viscous 
and  high  angle  of  attack  effects  and,  secondly,  the  propagation  of  noise  from 
the  rotor  blade,  particularly  the  nonlinear  propagation.  The  aerodynamics  work 
included  the  development  and  testing  of  a  Navier-Stokes  computational  solver 
for  rotor  blades  which  incorporates  rotating,  translating,  flapping  and  feathering 
motions.  Results,  which  focus  on  the  British  Experimental  Rotor  Programme 
(BEIRI*)  blade,  clearly  show  the  importance  of  including  all  motions  in  the  calcu¬ 
lation  of  aerodynamic  forces.  The  acoustics  research  concentrates  on  the  devel¬ 
opment  of  a  method  for  computing  the  nonlinear  propagation  of  acoustic  signals 
in  the  atmosphere.  The  method  is  based  on  a  boundary-element  discretization 
oi  the  time-dependent,  nonlinear  wave  equation.  Results,  computed  for  a  spheri¬ 
cally  symmetric  domain,  show  an  equivalence  with  VVhitham’s  method  up  to  the 
formation  of  a  shock. 
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1  Description  of  Research 

1.1  Aerodynamics  of  Advanced  Rotor  Blades 

This  section  summarizes  the  developed  Navier  Stokes  solution  method  which,  apparently, 
closely  resembles  that  employed  by  Duquc[l].  The  text  also  describes  the  grid  and  its  gen¬ 
eration.  The  current  grid  provides  a  notable  improvement  over  the  more  typical  C-11  grids 
utilized  in  previous  studies  of  the  BERP  blade,  since  it  can  more  completely  describe  the 
highly-swept  BERP  blade  at  the  tip,  a.  critical  part  of  the  geometry. 

Though  the  method  does  provide  improvements,  it  should  be  noted  that  some  incomplete¬ 
ness  remains.  The  first  problem  area  involves  the  wake  effect.  It  has  become  commonplace 
to  account  for  the  effect  of  the  rotor  wake  by  introducing  a  downwash  velocity  at  the  blade 
surface  as,  essentially,  a  boundary  condition.  This  surface  velocity,  computed  independently 
u.sing  a  code  such  as  CAMRAD,  can  be  uniform,  vary  radially,  or  vary  in  both  radial  and 
chordwi.se  directions,  depending  on  the  degree  of  complexity  sought.  However,  the  ])rograms 
which  generate  the  induced  velocities  cannot  at  present  handle  odd  geometries  or  high  blade 
angles,  and  are  of  little  use  for  predicting  the  induced  velocities  at  the  BERP  surface.  The 
current  method  ignores  the  effect  of  tlie  wake,  except  the  small  part  of  it  containe^l  within 
the  computational  grid.  It  is,  of  course,  impossible  to  predict  the  pressures,  forces,  and  mo¬ 
ments  which  would  arise  on  an  act  ual  lifting  blade  under  the  prescribed  flow  conditions  wlu'u 
the  wake  is  not  considered.  However,  there  appears  to  be  no  computationally  realistic  way 
to  include  any  reasonably  accurate  wake  effect  at  this  time.  Ignoring  the  wake  shouhl  iiot 
significantly  change  the  basic  characteristics  of  the  flow  field,  as  observed  by  'I'sung  r/.  ol.[2] 
who  use  only  a  constant  correction  factor  for  the  velocity  at  t  he  blade  surface.  ( !onsiderat  ion 
of  this  wake  effect  was  not  at  issue  for  i)u<|ue  who  studied  only  rectilinear  flight  regimes. 

File  remaining  area  of  concern  involves  the  turbulence  modeling.  It  appc'ars,  through 
a.  study  of  the  literature  and  through  discussions  with  ex|)erts  in  the  field,  that  tlu'rc'  rc'- 
mains  a  fair  amount  of  disagreement  regarding  the  accuracy  and  applicability  of  the  various 
turbulence  rnoch'ls.  It  seems  that  no  model  currently  exists  which  will  rc'liably  |)redict  flow 


Thin  Airfoil 


Cambered  Airfoil 


Reflexed  Airfoil 


Figure  1:  Rotor  airfoils  used  iti  the  current  study. 

separation  or  shock-wave  location  under  a  wide  range  of  conditions.  It  also  appears  in  doubt 
as  to  which  model  performs  best  under  the  conditions  under  study  here.  Based  on  some 
previous  results  of  Lomax  and  Mehta[.3]  and  Degani  and  Schiff[4],  and  because  of  its  relative 
simplicity,  the  Baldwin -Lomax  turbulence  model  was  chosen  for  use  in  this  study.  The  pri¬ 
mary  drawback  of  this  model  seems  to  be  in  predicting  shock  location  on  transonic  airfoils. 
The  previously  mentioned  references  indicate  that,  at  least  in  some  cases,  the  model  works 
well  in  predicting  separation  on  airfoils  and  bodies  at  high  incidence. 

d’he  advanced-geometry  rotor  blade  under  current  study  represents  only  an  approximate 
rendition  of  the  true  BERP  blade.  'Phough  the  blade  planform  is  available  in  the  puidic  do¬ 
main,  the  blade  airfoils  and  the  twist  of  the  rotor  inboard  of  the  tip  paddle  were  unavailable. 
I'he  rotor  studied  here  had  the  «;xact  BERP  planform,  luit  airfoils  were  defined  ba»efi  only 
on  the  sketches  i)ublished  by  Perry[.5].  Though  the  actual  BERP  blade  has  inboard  twist,  the 
current  blade  has  none.  Thus,  the  results  published  here  cannot  duplicate  the  performance 
of  the  actual  BERP  blade.  However,  since  one  objective  is  to  determine  planform  efh'cts  on 


the  flow  field  surrounding  the  rotor  blade,  the  actual  airfoils  used  should  have  only  secondary 
effect  on  the  desired  results.  Figure  1  shows  the  geometry  of  the  airfoils  used  in  generating 
the  results  presented  here.  One  set  of  calculations  was  performed  on  a  rotor  with  a  13ERP  tip 
but  with  NACA  OOxx  airfoils,  where  xx  refers  to  the  thickness  ratio  at  a.  particular  location. 
The  thickness  distributions  for  both  sets  of  airfoils  was  the  same. 

file  simplest  way  to  create  a  grid  for  solving  rotor-blade  aerodynamics  problems  consists 
of  generating  two-dimensional  C-grids  at  various  radial  locations  and  then  stacking  them 
in  the  spanwise  direction.  Imr  most  conventional  rotor  geometries,  this  approach  presents 
no  difficult}'.  Because  of  its  highly-tapered  tip,  howev('r,  using  this  approach  for  the  BFdlP 
blade  leads  to  extremely  small  cell  sizes  near  the  blad  tip  wliere  the  chord  l)ecomes  very 
smalt.  In  fact,  the  length  of  the  chord  on  the  act  u.ade  planform  goes  to  zero  at  the  tip — 
a  geometry  which  cannot  be  handled  by  the  C-H  generation  method  described  above.  An 
aiiproacli  to  avoiding  the  difliculties  at  the  tip  involves  “clipping"  the  blade  at  some  inboard, 
finite-chord  station.  This  approximation  can.ses  two  problems.  Firstly,  it  causes  the  loss  of 
the  effects  of  the  tip  regions  of  the  blade  which  may  be  of  great  importance.  ,Secondly,  even 
though  the  blade  is  clipiied  at  some  non-zero  chord  position,  for  the  geometry  to  Ik'  even 
closely  approximated,  the  chord-length  must  necessarily  be  (piite  small.  This  cause's  a  very 
(h'lise  clustering  of  grid  points  near  the  tip,  extending  radially  to  tlie  far  fiehl.  In  effect  ,  such 
a  grid  w,astes  grid  points  in  areas  where  they  are  not  needed  and  puts  a  severe  restrict  ion  on 
the  maximum  allowable  time  step  because  of  the  small  cell  size. 

The  problems  due  to  clipping  can  be  remedied  by  constructing  a  grid  in  which  t  he  “span- 
wis('”  stations  on  the  blade  surface  can  be  rotated  in  any  direct, ion,  specified  by  giving  leading 
and  trailing  edge  positions  of  a  station  rather  than  a  single  spanwise  location.  1  lu'  last  sta¬ 
tion  on  the  blade  can  thus  l)e  mad<'  parallel  to  tlie  actual  blade  tip.  Near  the  ti|).  both 
on  and  off  the  blade,  the  sections  are  specified  to  “fan”  into  the  inner  and  outer  parallel 
sections,  'fhe  upper  and  lower  surfaces  at  any  station  on  the  blade  are  interpolated  from  the 
actual  airfoils  using  biiiiu'ar  surface  patches,  hdgure  2  shows  a  computational  grid  on  the' 
blarle  surface. 

b 


Figure  2:  Surface  grid  on  the  BERP  blade  tip. 


The  only  equations  suitable  for  highly  separated  flows  and  amenable  to  presc'iit  com¬ 
puters  are  the  Reynolds-Averaged  Navicr-Stokes  eejuations.  An  implicit,  lime  accurate  LU 
decomposition  method  proposed  by  Obayashi  and  Fujii[6]  is  adapted  to  solve  for  the  flow 
properties  near  and  on  the  rotor  blade.  Benefit.s  from  this  procedure  include  less  computa¬ 
tional  work  as  well  as  less  required  storage  than  conventional  ADI  methods.  Furthermore,  it 
does  not  require  special  differencing  because  of  the  three-dimensionality  of  the  flow,  and  the 
code  to  implement  it  is  easily  vectorized.  A  local  geometrical  time  step  is  used  for  stf'ady- 
slate  problems  to  increase  the  convergence  rate.  Choosing  a  constant  time  step  will  produce; 
a  time-accurate  solution  for  unsteady  problems. 

rhe  solution  method  handles  the  blade  motion  by  imparting  a  vedocity  relative'  tee  a 
fixeel,  inertial  reference  frame  to  each  point  on  the  grid.  Grid  velocities  appear  in  the 
ge)ve*rning  equation,  df'rivc'el  in  a  fixed  frame  of  reference,  as  control  ve)lume  ve'locit ie's.  Feer 
heever  and  simj)le  forwarel  flight  cases,  the  expression  for  the'se  is  commonly  kne)wn  anel  fairly 
straightforward.  If,  however,  the  blade  undergoes  flapping  and  feathering  motions,  as  wc'll 
as  rotation  and  forward  flight,  the  transformation  from  a  grid-fixed  system  to  a  spacc'-fixed 
OIK'  becomes  much  more  complex.  Appendix  A  contains  a  derivation  of  the  grid  velocities 
relative  to  a  fixed  coordinate  system  for  a  pitching  and  flapping  rotor  blade.  .\t  this  time. 
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Experiment  □ 
FPR 

Viscous  Flow  - 


(a)v  =  30-  (b)v  =  60-  (c)v  =  90- 


I'igurc  3;  Coni[)a.risoti  of  curiviit  rosulfs  with  previous,  iiiviscirl  cah'ulatioti  aiirl  with 
('xix'rinu'tit  -  roctnngular  blade,  A//,,,  =  .7,  //.  =  .3. 

t  he  motion  is  constrained  to  tliat  for  a  rigid  rotor  [litching  about  any  radial  axis,  but  llap|)ing 
about  a  hinge  with  no  offset.  Including  a  hinge  offset  is  not  diflicult;  inclnding  (lexibility  in 
the  blade  requires  much  information  about  the  structural  modes  of  the  rotor  and  introdner's 
a  complexity  well  beyond  the  scope  of  the  current  si  iidv. 

1.1.1  Results 

Several  code  verification  ca.x's  were  run;  tho.se  shown  in  figure  3  compare  the  results  of 
t  he  current  code  with  thos('  computed  by  I''PH[7],  a  transonic,  full-pol('ut  ial  rotor  analysis 
code.  Some  ex[)erimental  data  is  also  inclurled.  I'he  figure  shows  pressnif'  ( oelficii'nl s  lor 
an  untwist('d  rectangular  rot.or  blade  with  a  constant  NA(.b\  0012  airfoil  at  /x'ro  angle  of 
attack,  with  rotational  tip  Mach  number  of  .7  and  advance  ratio  ..3.  'I'hree  azimut  hal  angles 
are  inebuh'd,  and  the  ladial  [)osition  is  at  89%  of  the  maximum,  d'lie  data  shows  good 
corrc'lation  wit  h  the  la'sults  of  bl’H  and  fairly  good  correlation  wit  h  tlu'  ('Xjx'rinu'iit . 

In  order  to  com[)ar<’  with  llu'  results  of  Diicpie,  the  code  was  run  for  the'  BI'.KP  planform 
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Original  computation 
Present  computation 
more  grid  points 


.0 

=  .2,  o  =  13°. 
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ill  uniform  flow  (non-rotating)  at  .2  Mach  number  and  13°  and  20°  angle  of  attack.  Figures 
4  and  5  show  the  computed  lift  coefficients  for  these  cases.  The  current  results  are  siiown 
for  two  grid  sizes.  The  original  grid  was  the  largest  possible  for  the  memory  available  on  the 
Arizuiia  State  University  Cray  X-MP.  When  the  NASA  Ames  Cray  Y-MP  became  available, 
a  grid  convergence  study  was  undertaken.  It  was  determined  that  an  increase  of  21%  in 
the  number  of  grid  points  was  required  for  convergence  with  the  blade  at  13°  angle,  and  an 
increase  of  81%  was  required  with  the  blade  at  20°.  The  maximum  grid  size  utilized  was 
123  x45  x43  points  in  the  wrap-around,  perpendicular,  and  spanwise  directions,  respect  ively. 

Several  discrepancies  exist  between  the  results  of  Duque  and  the  current  ones.  Uhese 
occur,  first  and  formost,  because  of  the  differing  airfoils.  Additionally,  since  Perry’s[r)]  and 
Duque’s  descriptions  of  the  airfoil  thickness  and  its  distribution  along  the  blade  differ,  some 
variations  in  the  lift  distributions  could  cc^ar  simiily  because  of  the  differing  thicknesses. 
This  would  be  particularly  apparent,  at  high  angles  of  attack.  Nevertheless,  the  major  trend 
of  lift  coefficient  with  span  is  consistent  between  the  two  results,  d’he  13°  case  was  run  with 
NACA  OO-T-T  symmetric  airfoils,  as  well  as  with  the  BF.RP-like  airfoils.  Results  from  these 
calculations  show  virtually  the  same  dependence  on  span,  but  at  a  slightly  lower  value  of  ry. 

A  major  thrust  of  the  work  dc'seribed  herein  is  to  determine  the  effect  on  the  advanced 
rotor  blade  of  the  unsteadiness,  not  only  due  to  forward  flight,  but  also  incurred  Ix’cause 
of  the  blade  motions.  Therefore,  the  Navier-Stoke.s  code  was  written  in  such  a  way  as  to 
easily  incorporate  arbitrary  blade  motion  through  time  dependent  grid  velocities,  .'\ppendix 
A  describes  the  derivation  for  the  grid  velocities  including  blade  feathering  and  flapping 
as  well  as  the  rotational  and  translational  motions.  'I'he  figures  b(4ow  show  results  for  a 
BFRP-like  blade  in  forward  flight  with  /i  =  .25  and  hover  tip  .Mach  number  equal  to  .65. 
'I’he  prescribed  collective  and  cyclic  i)itch  was  givim  as  =  10  -  lOsin  i/\  where  0  is  given  in 
degrees.  For  given  blade  properties,  it  is  (juite  possible  to  solve  for  the  blade  position  i.< 
the  blade  flapping  —as  a  conse(|uence  of  the  computed  forces  if  it  is  assumc'd  that  most,  of  tlu' 
lift  on  the  blade  comes  from  tlu'  outer  jnntion.s  over  which  the  comj)utation  is  carried  out. 
'I’liis  does,  however,  require  t  he  numerical  solution  of  yet  anot  her  diffi'i  iml  ial  ('quat  ion  which 
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describes  the  time-dependent  blade  motion.  In  order  to  avoid  the  additional  complexity 
(and  additional  computation  time)  of  including  the  blade  dynamics  at  this  initial  stage  of 
development,  the  flapping  angle  was  prescribed  to  be  (3  =  5  — 5  cost/’.  This  flapping  schedule 
is  not  necessarily  realistic,  especially  for  the  large  cyclic  pitch  introduced  above.  Nevertheless, 
it  does  provide  a  vehicle  for  examining  the  flow  response  to  the  type  of  unsteady  motions 
experienced  by  a  rotor  blade.  The  flapping  and  feathering  are  given  with  respect  to  the  shaft 
axis. 

Though  it  is  not  necessarily  valid  to  compare  the  results  from  the  rotor  in  forward  flight 
with  no  motion  to  those  from  the  rotor  with  flapping  and  feathering,  the  comparison  will, 
nonetheless  be  made.  As  shown  in  the  figures,  the  prescribed  motions  cause  large  changes  in 
the  pressure  distributions  on  the  rotor  blade  at  all  radial  stations.  The  illustrated  pressure 
distributions  are  shown  for  radial  positions  at  the  notch  and  in  the  middle  of  the  paddle. 

1.1.2  Conclusions 

The  developed  code  provides  a  robust  method  for  determining  the  forces  on  a  rotor  blade 
undergoing  general  motion.  The  current  code  could  be  enhanced  by  developing  a  manner  in 
which  to  include  the  effects  induced  by  the  wake,  a  task  which  cannot  be  completed  utilizing 
current  methodology.  Though  uncertainty  is  introduced  through  use  of  a  turbulence  model, 
little  can  be  done  to  ameliorate  this  in  the  absence  of  a  more  exact  closure  procedure. 

d’he  pre.sented  results  concur  with  some  previously  published  ones  and  introduce  some 
additional  conclusions  with  regard  to  the  aerodynamic  characteristics  of  the  BEHP  rotor 
blade: 

1.  At  high  angles  of  attack,  the  blade  |)roduces  a.  large,  stable  attached  vortex  at  the  ti]). 

2.  'I’he  notch  does  seem  to  play  a  role  in  BERP  aerodytiamics — a  vortex  is  produced  at 
the  notch  which  trails  into  the  wake.  This  vortex  does  not  appear  to  remain  attached 
to  the  blade,  although,  at  a.  20°  angle  of  attack,  there  may  be  soiik'  small  region  of 
attachment.  'I’he  notch  seems  to  ])rovide  a  stall  barrier  h*r  the  inboard  t  railing-c'dge 
stall  at  n  =  10°  and  for  t,h<'  large-scale  .separation  at  a  =  20°. 
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i5.  The  major  features  of  the  aerodynamics  of  the  paddle  section  are  unaffected  by  airfoil 
section. 

■t.  Feathering  and  flapping  motions  greatly  affect  the  loading  on  the  blade. 

1.2  Nonlinear  Propagation  of  Acoustic  Waves 

\Vhitham[8]  aiul  others  have  shown  t.hat,  in  situations  where  the  forces  due  to  acoustic 
pressure  disturbances  overshadow  those  caused  by  dissipative  ])roce.sses,  any  acoustic  signal 
will  experience  a  deformation  from  the  initial  waveform.  In  the  absence  of  dissipation,  the 
signal  will  eventually  steepen  and  form  a.  shock  due  to  the  nonlinear  nature  of  the  real 
atmosphere  in  which  the  signal  travels.  1  hough  tin*  amplitude  of  a  traveling  acoustic  wa\(' 
decreases  with  distance  from  the  source  of  the  disturbance,  the  deformation  of  thf'  wave 
becomes  more  pronounced  as  it  travels  farther  and  farther  from  the  position  of  its  inc('ption. 
Conseciucntly,  as  long  as  the  disturbance  is  loud  enough  at  the  source  (such  as  that  i)roduced 
by  a  fast-moving  rotor  blade),  the  signal  will  exhibit  nonlinear  steepening,  espc'cialh’  as  the 
wave  travels  to  the  far  field. 

Several  methods  arc  available  for  computing  the  wave  field  produced  by  an  acoustic  dis¬ 
turbance.  The  most  widely  atteinpti'd  involves  a  finite  difference  approach,  and  Ha('der[!)] 
has  enjoyed  some  success  using  this  method.  Finite  differencing  requiix-s  the  use  of  large, 
dense  grids,  and  it  tends  to  subdue  propagating  disturbances  because  of  the  inherent,  art  i¬ 
ficial  damj)ing  resulting  from  discretizing  the  partial  derivatives  of  the  governing  ('(piation. 
'Fhere  are  techniques  available  for  reducing  or  eliminating  this  extraneous  damping,  but  the 
current  research  centers  on  develoj)ing  a  methodology  for  solving  the  integral  ratln’r  than  tlu' 
differential  form  of  the  nonlinear  wave  equation.  Numerical  solution  of  tlu'  integral  ('(|uation 
involves  summations  rather  than  differences  ainl  is.  therefore,  less  subject  to  disert't izat ion 
errors  and  should  require  fewer  grid  points  for  an  accurate  solution. 

1  he  complete  acoustic  analogy  equation,  including  non-linear  effects,  can  be  generically 
written  as 

[./■(x,y,p(y,t))],_^  dF(y)  =  I),  (I) 


bl 


wIktc  p  roprc'sents  the  acoustic  (disturhanre)  density.  Strictly  speaking,  the  function  /  in 
the  integrand  depends  not  only  on  p,  hut  also  on  the  ’'ariahles  p  (acoustic  pressur(')  and  v 
(disturbance  velocity).  Additional  auxiliary  equations,  the  choice  of  which  d('i)ends  on  the 
assiiinptions  and  siinplincations  of  the  problem,  provide  the  relationships  betw('en  />.  v,  and 
p.  For  example,  the  isentropic  relation  provides  p  as  a  function  of  p  when  t  he  isentropic 
assumption  is  invoked.  The  otlu'r  governing  eciuations  give  velocity  for  one-dimensional 
propagation,  the  continuity  eciuation  proves  snfficic'iit  for  finding  tlu'  velocity  with  (huisity 
given. 

Equation  1  appears  similar  to  a  Volterra  integral  ecjuation  of  the  second  kind  but  diflers  in 
that  the  variables  in  the  integrand  must  be  evaluated  at  the  rvtnrdrd  t  ime,  r  =  /  —  |x  —  yl/oi, 
where,  as  is  customary,  x  repres('nts  ob.s<'rver  position,  y  represents  sourc('  position  (in  the 
integrand),  and  Tq  is  the  ambient  s|)eed  of  sound.  An  approximat  ion  to  eep  1  for  t  he  density 
at  N  ^-locations  (grid  si?/',  A'),  can  be  written  as 


-Y.d{pj{T))  [  /\'(x,y)dF(y)  =  0. 

j=i 


where  the  integrand  of  eep  1  is  replaced  with  //(p(y,  r))  A'(x,  y ),  and  i  runs  from  1  to  Ah 
I'hpiation  2  appears  very  similar  to  a  boundary  element  approach  to  solution  of  int('gral 
e([uations  except  that  the  ch'iisit.y  within  the  summation  must  Ix'  computed  at  the  r<-tardeti 
time.  Equation  2,  then,  represents  a  set  of  N  e<niations  (one  for  each  of  the  A’  grid  i)oints) 
with  N'^  unknowns — the  N  p,'s  and  the  i\'{N  —  I)  p^'s. 

In  order  to  make  the  comi)utational  procedure  more  tractabk',  a  simpler  form  of  e(j.  1 
was  utilized  to  solve  for  the  wa.v<'  fi<'ld  jiroduced  Ipv  a  pulsating  sphere: 
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Ec|uation  -d  is  efiuivalent  to  the  transonic  small-disturbance  ec|uation  for  aerodynamics,  and 
all  (luantities  have  been  made  dimensionless.  The  scpian'  bracl«'ts  indical  (' ('vahiat  ion  at  the 
retarded  l  ime.  In  ecj.  d,  g."  represents  the  velocity  |)<)t('nl  ial,  r„  is  the  normal  vclo<  ily  of  a 


l.fi 


surface  bounding  the  field,  and  II  is  tlie  distance  between  an  observer  point  and  a  point  in 


the  field.  'I'he  eejuation  can  be  discretized  for  an  N  x  ;U-point  grid  as 
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I'hpiation  5  represents,  effectively,  iV  x  M  e(juations  with  {N  x  M)^  unknowns,  just  as 
described  above  with  regard  to  eq.  2.  In  the  above,  n,j  is  the  angle  betwee-n  the  normal  to 
the  bounding  surface  and  the  vector,  x  —  y. 

droenenboomjlO]  s\iggests  using  a  linear  interpolation  of  the  unknown  varial)l(’  (in  this 
case,  9,)  as  a  function  of  /  in  order  to  estimate  its  values  at  the  appropriate  retarded  times. 
However,  the  linear  interpolation  for  9  gives  a  constant  value  for  d^/dl  for  the  durat  ion  of  a 
given  time  step,  At,  and  therefore,  gives  a  zero  value  for  the  second  derivative  with  respect  to 
t  except  at  times  given  by  iiAt  {11  an  integer)  when  tptt  is  undefined.  For  the  linear  probhun. 
in  which  (ftt  is  not  required,  this  simple  interpolation  may  work  quite  wc'll.  but  it  is  obviously 
not  suited  to  the  nonlinear  case.  It  was  decided,  therefore,  to  use  a  fast  h'ourier  transform  to 
interpolate  [)etween  known  times  to  determine  the  values  of  9  at  the  correct  retarded  t  imes. 
'I'he  known,  linear  solution  is  used  as  a  starting  .solution  for  the  iteration. 

Figure  8  shows  results  of  the  initial  computations  for  the  pressure  as  a  function  of  radial 
distance  for  a  Mach  number  on  the  surface  of  the  s[)here  of  .2  for  both  10  and  fiO  grid  points 
in  the  radial  direction.  I'he  figure  indicates  little  difference  bctwe('n  the  results  for  the  two 
grid  sizes.  Additionally,  the  figure  compares  tlie  results  to  the  solution  calculats'd  using  the 
classical  method  of  VVhitham  adaj>ted  to  the  current  problem  involving  a  continuous  wave 
train  in  three  dimensions.  I'he  computed  results  a|)|)ear  to  exhibit  the  classical  problem  with 
using  too  few  Fourier  coefficients  to  re.solve  (he  waveform.  An  attempt  to  solve  this  probh'in 
by  simply  (loubling  the  number  of  Fourier  coefficients  results  in  a  divergence  of  (he  iteration 
scheme. 

Using  the  Fourier  transform  (o  approximate  the  interim  .solution  after  each  iteration  has 
the  advantage  that  the  time  derivatives  in  eq.  fi  are  as  easily  estimat('d  ns  the  pottmtial 

If) 
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Figure  8:  Solution  at  t  =  0  for  four  Fourier  terms. 


function  itself,  and  no  numerical  approximation  to  derivatives  is  recpiired.  However,  each 
componeni  of  the  fburier  series  is  multiplied  by  its  frequency  when  the  first  derivative  is 
taken  and  by  the  frefpiency  squared  for  the  second  derivative.  I'his  procedure  tc'iids  to 
exagerate  any  errors  in  the  coefficients  of  the  higher  frequencies,  especially  in  the  nonlinear 
term  which  contains  the  product  of  the  first  and  .second  derivatives  of  the  velocity  i)otential. 
Consequently,  including  additional,  higher  frequencies  in  the  Fourier  approximation  to  the 
velocity  potential  leads  to  the  dominance  of  errors  in  these  frequencies.  As  the  nonlinear 
component  of  the  potential  becomes  more  j)ronounced,  such  as  when  the  Mach  number  of 
the  initial  disturbance  increases  or  the  distance  from  the  source  surface  becomes  large,  the 
importance  of  the  high  frequencies  and,  thus  of  their  errors,  also  increases. 

Figure  0  shows  the  computational  restdts  from  the  same  case  as  before,  but  with  eight 
terms  in  the  Fourier  series  and  using  a  technique  to  smooth  the  high  fix'fUK'iicie.s.  'flu' 
smoothing  method,  due  to  Laiiczos[l  1  ],  averages  the  value  of  the  signal  over  the  smallest 
miitTval  (A<),  or  highest  frefiuency,  in  the  series.  'Phe  technique  appears  to  damp  out  high- 
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Spatial  Pressure  Variation,  T  =  0 
M  =  .2,  8  Fourier  Components 


Figure  9:  Solution  at  <  =  0  for  eight  Fourier  terms  with  Lanczos  smoothing. 

freciueucy  errors  up  to  tlie  point  where  a  shock  (or  a  near-shock)  a|)pears  in  the  waveform. 
As  shown  in  the  figure,  the  computed  values  compare  well  with  the  Whitham  solution. 

1.2.1  Conclusions 

The  developed  computational  method  appears  to  provide  a  means  for  accurately  predicting 
the  nonlinear  propagation  of  an  acoustic  signal.  Because  the  nonlinear  effect  of  steepening 
becomes  progressively  more  ap|)arent  as  the  wave  travels  outward  from  the  source,  the  effect 
will  have  increased  importance  in  the  far  field.  Though,  as  the  figures  show,  the  method  gives 
good  results  for  .some  cases,  improvements  must  be  made  in  order  to  handle  more  general 
conditions.  Some  of  these  include; 

1.  Once  a  shock  (or  at  least  a  very  large  pressure  gradient)  appears  in  the  outward¬ 
travelling  wave,  the  Lanezos  smoothing  is  not  sufficient  to  stipress  the  error  in  the 
high  frequencies  of  the  Fourier  series  approximation  to  the  waveform,  d'herefore,  some 
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other  ineaiis  for  ('itlier  smoothing  the  Fourier  representation  or  even  for  a[)[)roxiniatitig 
the  function  must  l)e  attenijited.  A  suggested  smoothing  tecluncpu'  is  tlie  lilt('ring 
of  the  waveform  resulting  from  calculations  at  a  given  iteration  before  using  it  to 
compute  values  at  the  next  itcuation.  It  may  b('  more  appropriate  to  utilize  a  dilferent 
interpolation  method  altogether  for  d('t ('running  the  values  of  the  |)otential  at  tlu' 
retarded  times,  although  performing  an  Fl'd’  is  very  fast  relative  to  polynomial  sj)liii(’ 
fitting,  for  ('xample. 

2.  The  method  has,  at  present,  oidy  bcx'ii  tc'sted  for  the  very  simple  case  of  the  ]ndsat- 
ing  sphere,  which  is  both  stationary  and  spherically  symmetrical.  1  he  code  must  be 
adapted  to  account  for  more  general  geometries  which  do  not  lu'cessarily  exhibit  sim- 
{)lifying  symmetries,  and  for  moving  sources  which  arc  of  interest  in  the  calculation  (>f 
helicopter  rotor  noise,  d’hough  neither  of  the.se  adaptations  modifu’s  the  basic  method, 
which  is  a  solution  to  e(|.  2.  the  new  geometry  and  retarded  planform  calculations  add 
considerable  complexity  to  the  actual  computer  code. 

d.  The  iteration  procedure,  sometimes  termed  funciional  Hrrntioii,  seems  to  b('  i)ainfully 
slow  in  its  convergence  characteristics.  The  more  standard  methods  for  improving  the 
converg(’nce,  such  as  N'ewton-llaphson  or  secant  methods,  are  not  practical  for  such 
a  large  system  of  e(iuations.  Nevertheh'.ss,  to  Tuake  the  solution  mt'thod  feasible,  par- 
ticidarly  for  more  complicated  problems  with  larger  integration  grids,  its  convergence 
rate  should  be  eidianced. 

2  Supported  Personnel  and  Publications 
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the  .\avier-Stok('s  code,  was  suiij^orted  {(art -I ime  during  the  first  year  of  the  grant.  Dr.  Abdv 
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The  research  supported  by  the  grant  has  thus  far  resulted  in  one  publication: 


Abdy,  G.  and  Wells,  V.,  ‘b\  Numerical  Stiidy  of  Advanced  Rotor  Blades,’’  in  Pro¬ 
ceedings  of  the  International  Technical  Specialists  Meeting  on  Rotorcraft  Acous¬ 
tics  and  Rotor  Fluid  Dijnmnics,  American  Helicopter  Society  and  The  Royal 
Aeronautical  Society,  October  1991. 


It  is  anticipated  that  the  research  performed  under  this  grant  will  form  the  basis  for  several 
forthcoming  publications. 
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A  Derivation  of  Grid  Velocities  for  a  Rotor  Undergo¬ 
ing  Feathering  and  Flapping  Motions 

d'he  derivation  considers  a  rotor  blade  undergoing  four  basic  types  of  motion.  These  include 
1)  forward  flight  at  uniform  velocity,  2)  rotation  about  a  hub  a.xis  with  constant  angular 
velocity.  3)  feathering  about  an  axis  parallel  to  the  blade  span,  and  4)  flapping  about  a 
hinge  at  the  hub.  Because  the  code  reejuires  knowledge  of  the  grid  velocity  with  respect 
to  a  fixed,  inertial  coordinate  system,  the  above  motions  must  be  described  as  such  but  in 
terms  of  coordiates  on  the  grid.  Therefore,  if  Vj  represents  the  position  of  a  grid  j)oint  with 
respect  to  the  fixed,  inertial  system,  and  if  r'"  represents  the  position  of  such  a  point  in  the 
grid-fixed  system,  an  expression  is  required  which  gives  ry  in  terms  of  r'"  and  time.  Then 
the  time  derivative  of  ry  will  give  the  grid  velocity.  It  should  be  noted  that  this  exercise  is 
the  opj)osite  of  that  which  is  commonly  i)crformed  in  a  kinematics  j)roblem,  for  example, 
where  velocities  and  accelerations  are  generally  written  with  respect  to  a  body-fixed  system. 

On  an  actual  helicopter  rotor,  the  flapping  motion  arises  as  a  consequence  of  unequal  lift 
forces  at  various  azimuthal  angles  due  to  both  unsymmetrical  inflow  and  cyclic  control  inputs. 
'I’o  det('rmin<;  the  amount  of  flapping  and  the  aero<iynamic  forces  at  any  azimuth,  solution  of 
a  syst(mi  coupling  the  blade  dynamics  with  the  aerodynamics  is  recjuired.  Such  a  procedure  is 
far  too  time  consuming  for  the  current  task,  and,  for  this  reason,  a  flapping  motion  consideretl 
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reasonable  for  the  conditions  under  study  is  prescribed.  In  general,  flapping  and  featl.ering 


are  described  in  terms  of  Fourier  series  in  0,  the  azimuth  angle,  as 


0  =  0^^-^  A\  cos  0  +  sin  0  +  /I2  cos  20  +  B2  -sin  20  +  . . . , 


(5) 


(i  —  /?()  +  Cl  cos  10  +  61  sin  0  +  cos  20  +  b-z  sin  20  +  .  .  .  , 


(C) 


where  0  represents  the  pitching  angle  (featliering)  and  j3  the  flapping  angle.  The  current 
work  restricts  0  and  B  to  their  first  harmonics,  though  the  derivation  of  the  grid  velocities 
itself  is  general  enough  to  Include  any  specification  of  the  {)itching  and  flapping  motion. 

Consider  the  {x”' ,  y'” ,  z'")  coordinates,  shown  in  figure  lOe,  to  be  fixed  to  the  grid  (in 
this  case,  to  the  blade,  also)  such  that  y"'  extends  along  the  spanwise,  or  radial,  direction 
and  x'"  runs  in  the  chordwise  direction,  positive  toward  the  trailing  edge.  The  blade,  and 
grid,  pitch  about  an  axis  parallel  to  the  ?/"'-axis,  designated  by  the  ])osition  vector,  Fy.  Note 
that  the  vector,  Fy,  is  then  the  same  in  the  {x"\y"\z"')  system  and  the  {x",  y'\z")  system, 
whose  y-coordinate  coincides  with  y"\  but  which  does  not  pitch  with  the  blade.  Therefore, 
the  triple-[)rimed  coordinates  pitch  about  the  axis,  Fy,  with  respect  to  the  double-primed 
coordinates,  with  pitching  rate  0.  The  relationship  between  a  point,  f",  and  a  point  given 
by  r'"  is  then 

f"  =  Fy  -f  [//(/:)]  (r'"  -  Fy),  (7) 

where  n{0)  is  a  rotation  iiiatrix  between  the  two  systems.  Now,  the  {x'\  y",z")  coordinates 
flap  with  respect  to  a  system  designated  by  {x',y',z'),  whose  .c-axis  coincides  with  x'\  as 
shown  in  figure  lOc.  The  relationshii)  between  f'  and  f"  is  then 


F 


wheie  B(/i)  represents  the  flapping  rotation  ruatrix.  In  a  similar  manner. 


(!)) 


wlnu'e,  again,  .S'(0)  transforms  from  a.  rotating  (and  translating)  system  to  one  which  is  just 
translating.  (See  figure  10a.)  Note  that  the  {x,y,z)  system  is  translating  with  respect  to 
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fixed  space  at  a  constant  velocity,  V  which  may  approach  at  an  angle,  o,  so  that  tin;  velocity 
of  th('  {.r,y,z)  coordinates  with  r(;spect  to  the  fixcxl  ones  will  be  V  =  — V'cosox  +  I’siiwgiy. 
In  general,  then 

r^  =  r  +  V<.  (10) 

'I’he  /  snbscripts  refer  to  the  space-fixed  coordinates,  'lb  determine  the  grid  velocities,  note 
that 

ry-r  +  V.  (11) 

Now, 

r  =  lSU')mH)\  {ro  +  -  r„)),  (12) 

SO  that, 

r  =  [5Zi]r,)  1  p7?//](r"'-ro).  (Id) 

'I'he  rotation  matrices  are  all  fnnciions  of  time  according  to  the  i)rescribed  feathering,  flap- 
ping,  and  rotation  motions  which,  np  to  this  point  in  the  derivation,  can  be  con.sid('rcd 
perfectly  general.  'I'he  code  utilized  here  considers  only  motion  given  by  the  first  harmonics 
of  ecpiatittns  5  and  fi,  and  rotation  given  by  V’  = 

'I'lu!  code  which  [)erforms  the  flow  calctdalions  r<;qnires  only  the  current  grid  positions, 
r  or  ry,  and  not  those  in  the  grid-fixed  frame,  'lb  re<ince  the  storage  reejnirement-  -that  is, 
.so  as  to  not  have  to  store  r'" — the  code  solves  eq.  12  for  r'"  in  terms  of  r  at  eadi  time  step 
and  then  snl)stitntes  into  eep  l.'l  to  get  the  grid  velocities  in  terms  of  r  and  To  only. 
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